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Abstract 

Flux is an important source of uncertainties for a reactor neutrino experiment. It is determined from thermal power 
measurements, reactor core simulation, and knowledge of neutrino spectra of fuel isotopes. Past reactor neutrino 
experiments have determined the flux to (2-3)% precision. Precision measurements of mixing angle fl^ by reactor 
neutrino experiments in the coming years will use near-far detector configurations. Most uncertainties from reactor 
will be canceled out. Understanding of the correlation of uncertainties is required for flo experiments. Precise 
determination of reactor neutrino flux will also improve the sensitivity of the non-proliferation monitoring and future 
reactor experiments. We will discuss the flux calculation and recent progresses. 

Keywords: 

neutrino flux, reactor, fission rate 



1. Introduction 

Reactor neutrino experiments have played a critical 
role in the history of neutrinos. Among them. Savannah 
River Experiment 1 1 1 by Reines and Cowan in 1956 ob- 
served the first neutrino. Chooz ||2l determined the most 
stringent upper limit of the last unknown neutrino mix- 
ing angle sin2 26li3 < 0.17 in 1998. KamLAND IS ob- 
served the first reactor neutrino disappearance in 2003. 

Flux is an important source of uncertainties for a re- 
actor neutrino experiment. It is determined from ther- 
mal power measurements, reactor core simulation, and 
knowledge of neutrino spectra of fuel isotopes. Past re- 
actor neutrino experiments have determined the flux to 
(2-3)% precision. In the coming years, three precision 
experiments on neutrino mixing angle 9\t, using reac- 
tor neutrinos, Daya Bay |4|, Double Chooz |5|, and 
RENO ||6l, will start operation. All these experiments 
use near-far detector configurations. Most uncertain- 
ties from the reactor will cancel out. The residual er- 
ror will range from 0.1% to 0.45%. However, corre- 
lation among reactor cores need better understanding 
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of the error sources. Recently there are increasing in- 
terests on non-proliferation monitoring |7| using ton- 
level neutrino detectors. A reactor neutrino experiment 
at an intermediate baseline ~60 km with a giant detec- 
tor [8 1 will have rich physics content. Precise determi- 
nation of neutrino flux will greatly improve the sensi- 
tivity of these experiments. Another category of reactor 
neutrino experiments is v-electron or v-nucleus scatter- 
ing experiments, such as TEXONO |91, MUNU IfTOl . 
GEMMA UJJ, etc. Normally they won't rely on pre- 
cise neutrino flux. 

In this note, we will review the calculation of the re- 
actor neutrino flux and recent progresses on the error 
analysis. 

2. Calculation of Reactor Neutrino Flux 

Most commercial reactors are Pressurized Water Re- 
actor (PWR) or Boiling Water Reactor (BWR). They 
are very similar in neutrino flux calculation. We will 
use PWR as examples in the following. The ^^^U 
enrichment in fresh fuel of a PWR is normally (3- 
4)%, and more than 95% is ^■'^U. Electron antineutrinos 
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are emitted from subsequent yS-decays of fission frag- 
ments. They are dominated by 4 isotopes, ^■'^U, ^■'^Pu, 
^^'Pu, and ^-"^U. Other isotopes contribute only at 0.1% 
level. One can calculate the neutrino energy spectrum of 
each isotope by summing all fission fragment j6-decay 
branches. There are a lot of efforts on such studies. 
However, the fission products are very complex. Due 
to lack of accurate nuclear data, such calculations carry 
large uncertainties at the 10% level. The most accurate 
neutrino spectra of the first 3 isotopes were determined 
at ILL 1 12 1 by measuring the /3 spectra of fissioning. The 
/? spectra are then converted to neutrino spectra, with an 
average uncertainty 1.9%. The ^"'^U spectrum are calcu- 
lated theoretically lfT3l . They are shown in Fig.[T| The 
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Figure 1: Energy spectra of reactor neutrinos. 



isotope concentration in fuel will evolve during reactor 
operation as ^^^U depletes and ^^^Pu and ^"^'Pu breed. 
The ^^^U concentration is relatively stable. Such evolu- 
tion can be obtained by core simulation. A typical iso- 
tope evolution as a function of operation time, in terms 
of fission rates of the reactor, is shown in Fig.|2] 

When we know the fission rates of each isotope from 
core simulation, and neutrino energy spectrum of each 
isotope, we can easily get the neutrino flux S(Ey) - 

fiS iiEy), where fi is the fission rate of isotope / and 
S i(Ey) is its neutrino spectrum. However, the fission 
rates are proportional to the thermal power of the core, 
which is fluctuating. It is unrealistic to repeat core sim- 
ulation to reflect the power fluctuation. Normally we 
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Figure 2: Isotope evolution of a typical PWR. 

scale the neutrino flux to the measured thermal power. 



Zi(fi/F)ei 



J](fi/F)Si(E,), 



(1) 



where Wth is the thermal power, e, is the energy release 
per fission for isotope /, and F is the sum of thus fi/F 
is the fission fraction of each isotope. Among the in- 
puts, the thermal power data is provided by the nuclear 
power plant. The uncertainty is generally estimated to 
be (0.6-0.7)% L2,14J- Fission fractions are obtained by 
core simulation as a function of bum-up. Burn-up is the 
amount of energy in Mega Watt Days (MWD) released 
per unit initial mass (ton) of Uranium (TU). The sim- 
ulated fission fraction carries ~5% uncertainties from 
statistics of hundreds of analyses for various codes and 
various reactors j2irT5l. The 5% fission fraction uncer- 
tainties corresponds to ~0.5% uncertainty in neutrino 
yield. Energy release per fission varies slightly for dif- 
ferent cores at different time due to neutron capture and 
non-equilibrium products. Average numbers lfT6l [TTl 
can be used, with uncertainties of (0.30-0.47)%. Alter- 
natively we can extract them from the core simulation 
to accurately reflect the core differences and burn-up ef- 
fects. 

Recently there are studies to include contributions 
from non-equilibrium isotopes in the core ifTSl [T9l as 
well as that from spent fuel which is temporarily stored 
adjacent to the core ||20| |2Tl. These are at sub-percent 
level and only contribute to the low energy region. Be- 
sides fission products, ^^^U(n,7')^^^U reaction also con- 
tributes to the neutrino yield. It is below inverse yS-decay 
threshold (1.8 MeV) but will contribute significantly to 
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low energy v-electron scattering experiments ||9l. 

3. Thermal Power 

The most accurate thermal power measurement is the 
Secondary Heat Balance method. Detailed description 
of this measurement can be found, for example, in ll22ll . 
This is an offline measurement, normally done weekly 
or monthly. The uncertainty is cited as 0.7% by Chooz 
and Palo Verde. Primary Heat Balance tests are online 
thermal power measurement. Normally it is calibrated 
to the Secondary Heat Balance measurement weekly. 
Daya Bay power plants control the difference of these 
two measurements to less than 0.1% of the full power. 
These data are good for neutrino flux analysis. To 0.1% 
level, it can be taken as the Secondary Heat Balance 
measurement. The power plants also monitor the ex- 
core neutron flux, which gives the nuclear power. This 
monitoring is online, for safety and reactor operation 
control. It is normally calibrated to the Primary Heat 
Balance measurement daily. This measurement is less 
accurate, controlled to be less than 1.5% of the full 
power by Daya Bay power plant. 

Recently there are a lot of studies on the power uncer- 
tainties and instrumentation improvements by the power 
plants 1 031 122]| . with the motivation of power uprates. 
The power measurements can be more accurate than 
what were cited in the past reactor neutrino experi- 
ments. The uncertainties of the Secondary Heat Bal- 
ance is dominated by the flow rate measurement. In the 
past, there are two kinds of widely used flow meters, 
venturi type flow meters and orifice plate flow meters. 
Venturi flow meters are used by most US and Japan re- 
actors. The uncertainty is often 1.4%. It can be as low 
as 0.7% if properly calibrated and maintained. But they 
suffer from fouling effects, which could grow as high as 
3% in a few years. To improve the measurement, ultra- 
sonic flow meters have begun to be in use in some US 
and Japan reactors. They have uncertainties 0.45% for 
Type I and 0.2% for Type II |15|. The Orifice plate flow 
meters are used by French reactors. They have no foul- 
ing effects. Typically they have an uncertainty of 0.72% 
and could be improved to 0.4% with laboratory tests. It 
should be noted that the uncertainties of above flow me- 
ters are at the 95% C.L. (confidence level), as defined in 
ISO-5167. Unless specified, the thermal power uncer- 
tainty given by the power plant is also at 95% C.L. 

An EDF (Electricite de France) N4 reactor with four 
parallel steam generators, which is the Chooz type, is 
analyzed in l22\. Main components of the uncertain- 
ties are shown in Table [T] It is dominated by the dis- 
charge coefficient, which is an empirical formula in the 



flow rate measurement and its uncertainty is specified in 
ISO 5167-1-2003. The final uncertainty of the thermal 
power is 0.40% at 95% C.L. In this evaluation, it is as- 
sumed that the discharge coefficients of the orifice plates 
in four coolant loops are independent, thus the final un- 
certainty is statistically reduced. If the discharge coeffi- 
cients are fully correlated for all four orifice plates, then 
there is no statistical reduction. The power uncertainty 
will be 0.37% at Icr level, which is still significantly 
smaller than 0.7% that Chooz used. 



Table 1 : Error table for the thermal power measurements for N4 reac- 
tor (22). 





Contribution 


Relative fraction 


Origin of uncertainty 


[MWth] 


[%] of the 






17.2 MWth 


Discharge coefficient 


15.33 


79.57 


Differential pressure 


6.33 


13.57 


Steam gen. inlet temp. 


2.81 


2.68 


Primary input 


2.00 


1.35 


Others uncertainties 


2.98 


3.00 


Uncertainty at 95% C.L. 


4250±17.2MW (0.40%) 



The Daya Bay and Ling Ao reactors are all calibrated 
with the SAPEC system, which is an EDF portable high 
precision secondary heat balance test system with its 
own sensors, databases, and data processing, of uncer- 
tainty of 0.45%. The calibration results of Ling Ao reac- 
tors can be found in Table[2][23J. Four tests show differ- 
ences from 0.031% to 0.065%. Such small differences 
mean that the secondary heat balance system (KME) 
of the reactor is strongly correlated with the SAPEC 
system. Actually, the calibrations use the same orifice 
plates but different pressure transmitters. It proves again 
that the power uncertainty is dominated by the discharge 
coefficient. It also shows that Ling Ao KME is in very 
good agreement with the SAPEC system. The power 
uncertainty is estimated to be 0.48% at 95% C.L. in this 
comparison. As with the above example, these analyses 
assumed that the discharge coefficients are uncorrelated, 
which may not be the case. 

4. Core Simulation 

The fission fraction of fuel isotopes are obtained by 
core simulation. Qualified core simulation codes are 
normally licensed, and not available to scientific col- 
laborations. The core simulation also needs a lot of in- 
formation from the power plant as inputs. Fortunately, 
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Test 1 


Test 2 


Test 3 


Test 4 




KME(MW) 


2897.1 


2904.4 


2908.9 


2906.9 


Thermal 


SAPEC(MW) 


2896 


2903 


2907 


2906 


Power 


Difference(MW) 


1.1 


1.4 


1.9 


0.9 




Difference 


0.038% 


0.048% 


0.065% 


0.031% 


Uncertainty 


KME 


0.4806% 


0.4806% 


0.4806% 


0.4806% 


Analysis 


SAPEC 


0.45% 


0.45% 


0.45% 


0.45% 



the fission fraction can be extracted as a by-product of 
the refueling calculation required by the power plant, 
as a function of bum-up. The uncertainties of the ob- 
tained fission fraction depends on the simulation code. 
It only slightly depends on the inputs such as tempera- 
ture, pressure. Boron concentration, etc. 1,24] , as tested 
by the simulation code ROCS. The uncertainties of the 
simulation can be studied by comparing the measured 
and calculated concentration of fuel isotopes sampled 
at different bum-up when refueling. These studies are 
normally a part of the qualification of the licensed sim- 
ulation code. In ref. Ll_5|, 159 such studies for vari- 
ous codes and various reactors in US and Japan haven 
been collected and analyzed. On average, the simulated 
concentration of isotopes have uncertainties of ~4% for 
235u, -5% for 239pu, -6% for 24ipu, and -0.1% for 
2^^U. Assuming the simulated neutron flux in the core 
is not affected by small variations of isotope concen- 
tration, the fission rate is proportional to the isotope 
concentration. Due to the strong constraint of the to- 
tal thermal power, simulated concentrations of the foure 
isotopes are not independent. A 5% error on the iso- 
tope concentration corresponds to a ~0.5% uncertainty 
on the detected neutrino rate via inverse y6-decay reac- 
tion. There also are large core-to-core correlations for 
reactors simulated with the same code. 

The energy release per fission e,- in Eq. 1 is defined 
as the energy absorbed in the reactor per fission event. 
Neutrinos will take away energy from the total energy 
released in a nuclear fission. Some fission fragments 
of long lifetime will not reach equilibrium in a short 
time, thus part of the energy is not released. Each fis- 
sion will produce 2-3 neutrons while only one neutron 
will be used to maintain the chain reaction at stable 
state in reactor running. Other neutrons are absorbed 
and release energy via neutron capture. As the fuel 
composition evolves, the contribution of neutron cap- 
ture will change. Taking the above three connections 
into account, e, varies slightly for different reactors and 
it changes with time. Average e,'s were evaluated for 



typical reactor configurations in refs. |16| and |17|. A 
more accurate estimation of e,- as a function of burn-up 
can be obtained from core simulation, where these cor- 
rections are included automatically. 

Table 3: Energy release per fission in MeV from refs. 1161 and 1171 . 



Isotopes 


James 


Kopeikin 


235U 


201.7±0.6 


201.92+0.46 


238u 


205.0±0.9 


205.52+0.96 




210.0±0.9 


209.99+0.60 


Pu 


212.4±1.0 


213.60+0.65 



5. Neutrino Spectra 

Electron antineutrinos are emitted from the subse- 
quent j6-decays of fission fragments. Due to the lack of 
data for the yS-decays of the complex fission products, 
theoretical calculations of the neutrino spectra of iso- 
topes carry large uncertainties. ILL [12J measured the 
/3 spectra of fissioning of ^^^U, ^^^Pu, and ^*^Pu by ther- 
mal neutrons, and converted them to neutrino spectra. 
The normalization error is estimated to be 1.9%. Spec- 
trum shape error is from 1.34% at 3 MeV to 9.2% at 8 
MeV, as shown in Fig.[r| ^^^U can not fission with ther- 
mal neutrons. Its spectrum relies on theoretical calcula- 
tion. The uncertainty is estimated to be 10% |13|. Nor- 
mally 238u contributes (7-10)% of fissions in a PWR. 
The calculated neutrino counting rate and spectra were 
verified by Bugey and Bugey-3 1251 . The normaliza- 
tion error is further lowered to 1.6%, which was used 
by Chooz. 

6. Non-equilibrium Isotopes and Spent Fuel 

The ILL spectra are derived after 1.5 days exposure 
time with thermal neutron. Thus, long-lived fission 
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fragments have not reached equihbrium. In a real re- 
actor, these fission products will accumulate and con- 
tribute to the neutrino flux. Chooz estimated this contri- 
bution to be ~0.3% on the average and ignored it in the 
detailed analysis due to its small size comparing to other 
errors. Six chains have been identified in L18J . with half 
lives from 10 hours to 28 years. They only contribute to 
the low energy region. Further studies show that for a 
typical PWR, on average these contributions are ~Q.2% 
of total neutrino detection rate via inverse yS-decay lfT9ll . 
In the 2-4 MeV region, it increases to 0.8% after one 
year's accumulation, as shown in Fig. [3] 
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Figure 3: Contribution of non-equilibrium isotopes as a fraction of the 
total neutrinos, weighted by the inverse /?-decay cross section tl9l . 



Spent fuel is normally stored temporarily adjacent to 
the core. The storage could be as long as 10 years. Sim- 
ilar to the non-equilibrium contributions, the long-lived 
fission fragments in the spent fuel will contribute to the 
neutrino flux. A PWR is normally refueled very 12-18 
months. The spent fuel from one refueling will con- 
tribute ~0.2% of the total neutrino rate after first several 
days. 

7. Conclusion 

Before the 1980's, the reactor neutrino flux was de- 
termined to an uncertainty of 10%. With a lot of ef- 
forts, especially by ILL, Bugey, Chooz, Palo Verde, 
etc., the uncertainty has improved to (2-3)%. Motivated 
by high precision neutrino measurements by scientific 
collaborations and power uprates by the power plants, 
we have more accurate thermal power The uncertainty 
could be lowered from 0.7% to 0.4%. Small correc- 
tions from non-equilibrium isotopes and spent fuel, as 
well as energy release per fission are studied in detail. 
We also have a global picture of uncertainties of fission 
rate simulations. However, there is no new data on neu- 
trino spectra of fuel isotopes. For single-detector exper- 
iments, the neutrino spectrum uncertainty of about 2% 
will dominate. The next 6*13 experiments with near-far 



relative measurements will suffer little from reactor flux 
uncertainties, which is estimated to be 0.1% to 0.45%, 
depending on the experiment layout. The correlation 
among reactor uncertainties is important for fl^ exper- 
iments since correlated errors will cancel out. Mean- 
while, high precision detector and high statistics at the 
near detector of these experiments may help to improve 
the knowledge of neutrino spectra. 

Since the time that this manuscript was written, it has 
been suggested ll26l that the antineutrino spectra for all 
of the relevant fission isotopes (^"'^U, ^"^^U, ^^^Pu, and 
^^'Pu) are several percent larger than those of ref. IIT2I 
[T3l . This issue requires further attention. 
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